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A new over-improved stout-link smearing algorithm, designed to stabilise instanton-like objects, 
is presented. A method for quantifying the selection of the over-improvement parameter, e, is 
demonstrated. The new smearing algorithm is compared with the original stout-link smearing, and 
Symanzik improved smearing through calculations of the topological charge and visualisations of 
the topological charge density. 
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I. INTRODUCTION 

Studies of long distance physics in Lattice QCD sim- 
ulations often require the suppression of short-range UV 
fluctuations. This is normally achieved through the ap- 
plication of a smoothing algorithm. The most common 
prescriptions are cooling [J, 0, [|[ , APE [JJp] , and im- 
proved APE smearing [y], HYP smearing Q and more 
recently, EXP or stout-link smearing [8| and LOG smear- 
ing [9| . Filtering methods such as these are also regularly 
used in calculations of physical observables to improve 
overlap with low energy states. 

All smoothing methods require the use of an approxi- 
mation to the continuum gluonic action 



s - 1 - 



d A x tr[F^F^} . 
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Because we are approximating space-time by a 4-D lat- 
tice, these approximations will contain unavoidable dis- 
cretisation errors. These errors can have a negative effect 
on the topological objects present in the gauge field being 
studied and are detrimental to the smoothing process. 

Successful attempts have been made to remove the low 
order discretisation errors by combining Wilson loops of 
different shapes H, fiol. [Til [~L2| . Unfortunately, as shown 
by Perez, et al. [LI and reiterated below, it is still possi- 
ble for the remaining errors to spoil instanton-like objects 
in the field. They proposed over-improved cooling as a 
means of taming these errors via the introduction of a 
new tunable parameter e into their action p~3| . 

In Sec. [IT] we begin by presenting a summary of the 
most common smoothing algorithms. We then motivate 
the introduction of over-improvement by considering the 
effect of the discretisation errors on the action of a sin- 
gle instanton in Sec. IIIII Extending the work of Perez, 
et al. [l| we define in Sec. lIVI an over-improved stout-link 
smearing algorithm and demonstrate how to quantita- 
tively select a value for the parameter, e. Lastly, in SecfVl 
with the calibration of the over-improved parameter com- 
plete we test the over-improved stout-link smearing on a 
variety of lattices, including large 28 3 x 96 light dynam- 
ical gauge fields from the MILC collaboration [Tl.TT5j. 



II. SMOOTHING ALGORITHMS 

Standard cooling proceeds via a systematic sequential 
update of all links U^(x) on the lattice, where at each link 
update the local Wilson action Hi is minimised. The 
local Wilson action corresponding to U^{x) is defined as 



S w {x) 
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where 



= U v (x + A) Ul(x + v) Ul{x) 
+ Ul(x + fj,-v) Ul(x - v) U u (x - 0) 



(3) 



is the sum of the two staples touching Un{x) which reside 
in the [i-v plane. From Eq. ([2|), we can see that Sw 
will be minimised when Retr [1 — U ^{x)Yi ^ v (xj\ = 0. It 
naturally follows that when cooling, the aim is to replace 
Ufj,{x) with a new link that optimises 



max Re tr 



U^x) 

V 



(4) 



When performing this update in parallel, one must be 
careful not to replace any link which is included in the 
local action of a neighbouring link. This requirement 
means that cooling is a relatively slow operation com- 
putationally, but fast in regard to the removal of action 
from the gauge field. 

APE smearing differs from standard cooling in that all 
links can be simultaneously updated in a single sweep 
through the lattice, resulting in a significant speed in- 
crease. In APE smearing one first calculates a smeared 
link U'(x), which is the weighted sum of its nearest neigh- 
bours. 



U'Jx) = (l-a)U li (x) + %YxUx), 



(5) 



where E Mt/ is defined as in Eq. ([3]) , and a is a real param- 
eter, usually set to ps 0.7. The new link U'{x) is then 
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projected back into the SU(3) group via some projection 
operator V, 

U„(x)=VU' fl (x). (6) 

The projection of Eq. (J6j) is necessary because we have 
performed an additive step in Eq. ([5]), which is not 
an SU (3) group operation. The projection step is not 
uniquely defined, but the preferred method is to select 
the new smeared link Uu(x) such that it maximises 

Retr(tf M (aOE#(aO) • (7) 
In the limit a — > 1 we see that Eq. ([5]) becomes 

W-^E 53 ]^)- (8) 

UjtfJ, 

Substituting this result into Eq. ([7]) shows how the pro- 
jection method has become equivalent to cooling Eq. Q, 
and that there exists a direct link between APE smearing 
and cooling in the limit that links are updated sequen- 
tially. The simultaneous update of APE smearing limits 
a < 0.75 [l7|. 

The more recent smearing technique, stout-link smear- 
ing Q , makes use of the exponential function to remain 
with the gauge group and remove the need for a projec- 
tion step. Beginning with the staples Eq. define 

where we will choose an isotropic four-dimensional con- 
stant p^u = p sm , but other selections are possible. The 
matrix defined by 

Q^) = ^l(x)-%(x)) 

-itr^x)-^)), (10) 

with 

%(x) = C^(x)Ul(x), (11) 



is by definition Hermitian and traceless, and hence 
g^O) g SU(3). The new smeared link is then defined 
by 

U fl {x)=exp(iQ ll (x))U li {x). (12) 

An expansion of the exponential in Eq. (I12j) results in 
the same sum of paths, to first order in p sm , as for APE 
smearing Q. Given this, and the already established 
link between APE smearing and cooling, it follows that 
there exists a connection between cooling and stout-link 
smearing. Indeed, simulations of Lattice QCD show that 
for any given gauge field, the structures revealed by the 
smoothing procedures are remarkably similar. 



III. DISCRETISATION ERRORS AND 
IMPROVEMENT 

The corrosion of topological structures during the 
smoothing process is a well known side-effect of both 
cooling and smearing [1, [HI, Q3 • It is the unavoidable dis- 
cretisation errors in the lattice action that are the cause 
of this observed behaviour. This obviously inhibits our 
ability to study the topological excitations on the lattice 
with the most local operators and it would be beneficial 
if it could be prevented. 

When a gauge field is smoothed, the topological struc- 
tures within are subjected to the effects of lattice dis- 
cretisation errors. One such topological excitation is the 
instanton. To understand how the errors will alter instan- 
ton distributions, first consider the clover Wilson action 
given by 

S W = 0^52 (l-P^x)) , (13) 

where P^ v (x) denotes 1/3 of the real trace of the clover 
average of the four plaquettes touching the point x. 

Following Ref. [lU, Sw can be expanded in powers of 
the lattice spacing, a, giving: 



Sw — fl4 J]] J]] tr 



X [1>V 



\F%,{x) - a - {{V,F, v {x)f + (V v F, v {x)f) 



1 



1 



24 

+ 0(« 10 ,/), 



(14) 



where igF^ = [D^A,], A = + igA^, and T> tl 4> = <j>], for arbitrary <j>. 
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The goal is to substitute the instanton solution [It 
given by 



IV. OVER-IMPROVEMENT 
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where 



S 



Xi ± i x ■ a 



(15) 



(16) 



for instantons/anti-instantons with a the Pauli matrices, 
into the expanded Wilson action Eq. (fTj|) . This requires 



the use of the lattice approximation a 4 ~ / d 4x - Sub- 



stituting the instanton solution Eq. (| 1 5[) into Eq. ([14]) and 
performing the integration then yields, 
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70 
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(17) 



Notice that the leading error term in Eq. (fTF)) is nega- 
tive and depends upon the instanton size Pmst- When 
the Wilson action is used in a smoothing algorithm these 
errors result in an under- estimation of the action density. 
Additionally, by decreasing pi ns t the action will be fur- 
ther reduced. The smoothing algorithms, which are try- 
ing to decrease S, will therefore shrink pi ns t in order to 
reduce the action. Repeated application of the smoothing 
procedures will eventually lead to overwhelming discreti- 
sation errors and cause instantons to "fall through the 
lattice" and disappear. 

Improved actions aim to fix the problem of discretisa- 
tion errors by including different sized Wilson loops in 
the calculation of the action, S. By choosing the coeffi- 
cients of the loop combinations carefully it is possible to 
eliminate the leading order error terms. 

The Symanzik improved action uses a linear combina- 
tion of plaquette and rectangular loops to eliminate the 
0(a 2 ) errors. 



-{l-P^x)) 



--((l- RfW (x)) + (l-B^(x))) 



(18) 



Analogous to , i? M „ and R vj]b denote the different pos- 
sible orientations of the rectangular loops. 

This can be expanded in terms of a, and the instan- 
ton solution substituted as above to find, in agreement 
with 111 , that 
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(19) 



The 0(a 2 ) error term has been removed by design, but 
we see that the 0(a 4 ) term is still negative. Therefore, 
this action will still not preserve instantons. 



Formalism 



In 1993, Perez, et al. [l3j introduced the notion of over- 
improved cooling, also known as e-cooling. The essen- 
tial idea was that instead of trying to use different loop 
combinations to completely eliminate higher order error 
terms, they would instead choose their coefficients such 
that the error terms become positive. 

Introducing the parameter e, they defined the following 
action, 



4-e 



+ 6 -^-{l-W, u {x)) 



(20) 



where W^(x) denotes the clover average of the 2x2 
squares (windows) touching the point x. Note that 
in Eq. (f2"0"l) , e has been introduced such that e = 1 gives 
the standard Wilson action and e = results in an 0(a 2 ) 
improved action. Expanding Eq. (I20|) in terms of a and 
substituting the instanton solution Eq. (fT5|) gives 
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(21) 

where the 0(a 2 ) term is directly proportional to — e. 
Thus, by making e < the leading order discretisation 
errors become positive, and the modified action should 
preserve instantons. 

In the interests of preserving locality we choose to use 
the traditional combination of plaquettes and rectangles 
as in the Symanzik improved action in preference to the 
combination of the lxl and 2x2 loops used in [l3| • We 
now introduce the parameter e into the Symanzik im- 
proved action Eq. (fT8"|) . By requiring that e = gives the 
0(a 2 ) improved Symanzik action, and that e = 1 gives 
the standard Wilson action. This implies the following 
form for the action, 



1 - e 



5 - 2e 



12 



((1 - R^{x)) + (1 - JM*))) 



(22) 



Performing the usual expansion in a gives: 
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(23) 



Into which we substitute the instanton solution to find 
that 
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As in Eq. (|21[) . negative values of e will result in a positive 
leading error term, and should preserve instantons. 

We introduce the over-improvement parameter into 
the stout-link smearing algorithm by modifying the link 
combinations used in Eq. ©. Whereas the original 
C n (x) = p sm XX 1 x 1 paths touching C/ M (x)}, the mod- 
ified stout-link Ca(x) has the form 



E 



2c 



(lxl paths touching U u (x)) 



12 



■(1x2 + 2x1 paths touching U^x)) 



(25) 

and the smearing parameter p sm is unchanged. Note 
that both forward and backward horizontally orientated 
rectangles are included in the 2x1 paths, such that il u (x) 
resembles the local action. 



B. Tuning 

Of course, this now begs the question: How negative 
should e be in order to preserve instantons? Perez, et 
al. reported a value of e = —1 to preserve instantons, 
and indeed it does. However, just as positive values 
of e can shrink instantons, so too can negative values 
cause instantons to grow. Just as small instantons can 
fall through the lattice, big instantons can grow so large 
that they are destroyed by the smoothing procedure [19j | . 
Additionally, one does not want to unnecessarily distort 
the instanton-like objects in the gauge field. Care must 
therefore be taken not to choose a value of e that is too 
negative. 

In order to quantify the selection of e we consider the 
ratio S(e)/So, where So = 87r 2 /g 2 is the single instanton 
action. Ideally S(e)/Sq should be equal to 1 for all values 
of the instanton size, pinst-, as it is in the continuum. 

Plots of S(e)/So versus pi nst for the Wilson and 
Symanzik actions are shown in Fig. [T] Note that it is ac- 
tually the slope of the curve that will govern whether an 
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FIG. 1: S(e)/So versus the instanton size for the Wilson and 
Symanzik improved actions. The ideal action would produce 
a flat line at S(e)/So = 1. The positive slope on both curves 
means that instantons will shrink if the Wilson or Symanzik 
actions are used to smooth the gauge field. 



instanton shrinks or grows, not just the sign of S(e)/So- 
Although the Symanzik action is closer to the ideal ac- 
tion than the standard Wilson action, the slope is still 
positive for all Pi ns t and using this action will shrink in- 
stantons. 

The goal is now to select a value of e that results in 
the flattest line possible, whilst ensuring the stability of 
instantons. A plot for three different values of e is shown 
in Fig.O With e = — 1 the curve is similar to the mirror 
image of the Wilson action. For api ns t > 1.5, e = —0.25 
and —0.35 give curves closer to the ideal, however as |e| 
is decreased the maximum occurs at larger pi ns t- Since 
it is the slope that is responsible for how an instanton's 
size changes, the maximum of S(e) gives the dislocation 
threshold of the smearing algorithm. Assuming that any 
topological excitation of length > 2a is not an unphysical 
UV fluctuation or a lattice artifact, one should aim for a 
dislocation threshold of < 2a. 

Given this, we propose that a value of e = —0.25 will 
be sufficient. This choice gives a dislocation threshold of 
~1.97, and a curve that is mostly flat down to values of 
Pinst ~ 1-7. The action S(e)/So is also very close to the 
ideal. 

In Fig. [3] we provide a comparison of the Perez over- 
improved action, our over-improved action S(— 0.25), and 
the standard Wilson action. It is clear that S(— 0.25) 
should produce the best results. 
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FIG. 2: S(e) for three different values of e. The larger — e is 
made the further the curve moves from the ideal behaviour 
and the sharper the maximum. 
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FIG. 3: A comparison of S(e)/So for the Perez over-improved 
action, our over-improved stout-link smearing S(— 0.25), and 
the standard Wilson action. 



Given a value for e one can find a suitable value for the 
smearing parameter, p sm . Starting from some arbitrary 
value, systematically increase p sm until uq (the mean- 
plaquette value) no longer increases when smearing. This 
value sets an upper threshold for p srn and one should then 
choose some p sm suitably below this threshold. In what 
follows we use a value of p sm = 0.06. A typical value 
for standard stout-link smearing is p sm « 0.1. The over- 
improved algorithm is more sensitive to the smearing pa- 
rameter than standard smearing because of the larger 
loops used in the smoothing procedure. 

We mentioned earlier that it was the slope of S (e) / So 
that is responsible for how an instanton will evolve under 
smearing. We thought it prudent to check that this was 
actually the case in practice by smearing a single instan- 
ton gauge field. In order to exaggerate the effects, we 
selected a rather extreme value of e = — 4. A comparison 
of e = —4 and e = —0.25 is given in Fig. [4]. Using e = —4 
to smear a single instanton of size pi ns t ~ 1 should de- 
stroy the instanton. Meanwhile, if we smear an instanton 
with size pi ns t ~ 3.5 it should grow rapidly when e = —4, 
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FIG. 4: The evolution of a single instanton under over- 
improved stout-link smearing with e = —4.0 and e = —0.25. 
The S(e)/So plots for the respective e values are shown in 
the top plot. The bottom plot shows how an instanton of 
size pinst ~ 3.5 grows much more rapidly with e = —4.0 
(upper solid line) than for e = —0.25 (dashed curve). The 
lower solid curve in the bottom plot is for an instanton of size 
Pinst ~ 1 and e = —4.0. It has been offset by +1.8 to fit 
on the same scale as the larger instanton. We see how the 
instanton initially shrinks and is destroyed after 20 sweeps of 
over-improved stout-link smearing with e = —4.0. 



and stay relatively the same size for e = —0.25. 

Fig.0]shows the results of the simulations. The lowest 
curve is for the small instanton. Note that we offset by 
the size of the smaller instanton by +1.8 in order to dis- 
play it on the same scale as the larger instanton. As ex- 
pected, the small instanton initially shrinks under smear- 
ing until it is eventually destroyed in the 20 th sweep. For 
the larger instanton we see that it does grow rapidly for 
e = —4.0, but that its size stays relatively constant for 
e = -0.25. 



V. ALGORITHM COMPARISONS 

Given the selection of e = —0.25 it is now important 
to make a comparison of over-improved stout-link smear- 
ing with normal stout-link smearing. We are primarily 
concerned with the stability of the topological charge un- 
der smearing, and the structure of the gluon fields after 
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smearing. 

We use two sets of gauge fields for this study. Firstly, 
an ensemble of large 28 x 96 dynamical MILC lattices [14|, 
ITsT ]. with light quark masses; am Ui d — 0.0062, am s — 
0.031. We will also use a quenched MILC ensemble of the 
same size and lattice spacing a = 0.09. The gauge fields 
were generated using a Tadpole and Symanzik improved 
gauge action with 1x1 + 1x2 + 1x1x1 terms and 
an Asqtad staggered dynamical fcrmionic action for the 
2 + 1 flavours of dynamical quarks. 

We also use some quenched CSSM gauge fields cre- 
ated with the 0(a 2 ) mean-field improved Liischer-Weisz 
plaquette plus rectangle gauge action [20| using the pla- 
quette measure for the mean link. The gauge-field pa- 
rameters are defined by 

5/3 



S G =^£(l-P M „(z)) 



£ (2-iM*)) 



(26) 
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The plaquette measure of the tadpole improvement factor 
is 



/ \ 1/4 

u = [(PnA%)) x<fl<v ) 



(27) 



where the angular brackets indicate averaging over space- 
time and plaquette orientations. The CSSM configura- 
tions are generated using the Cabibbo-Marinari pseudo- 
heat-bath algorithm [2l[ using a p arallel algorithm with 
appropriate link partitioning [221 ] - To improve the er- 
godicity of the Markov chain process, the three diagonal 
SU(2) subgroups of SU(3) are looped over twice [6| and 
a parity transformation [231 ] is applied randomly to each 
gauge field configuration saved during the Markov chain 
process. 
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A. Topological Charge 
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Let us first consider the evolution of the total topo- 
logical charge of a gauge field under stout-link smearing. 
Typical studies in the past have rated a smearing algo- 
rithm's success by its ability to generate and maintain 
an integer charge. We will also use this test to rate the 
effectiveness of the smearing procedures because of its 
simplicity and widespread use. However, it should be 
noted that we will be smoothing extremely large 28 3 x 96 
lattices. Due to the vast amount of non-trivial topolog- 
ical charge field fluctuations present it will take a lot of 
smoothing to generate an integer charge. 

Fig. [5] provides a sample of 4 different gauge fields 
smeared by standard, Symanzik improved, and over- 
improved stout-link smearing. The first is a 28 3 x 96 
quenched MILC gauge field, the centre two are 28 3 x 96 
light dynamical MILC fields, and the last is a smaller 
16 3 x 32 quenched field. 



FIG. 5: Plots showing how the topological charge evolves 
under stout link smearing when using three different actions. 
The top graph is for a 28 x 96 quenched gauge field. The centre 
two are from an ensemble of 28 x 96 dynamical fields with light 
quark masses. The bottom is a smaller 16 x 32 quenched gauge 
field. The features of the graphs are explained in the main 
text. 



The top graph shows an example of the over-improved 
action producing a stable result. In this instance the Wil- 
son action is fluctuating widely, and is unable to reach 
a stable charge within 200 sweeps of smearing. The 
Symanzik improved action is better in that it stabilises at 
around 120 sweeps, however the over-improved action is 
clearly superior, stabilising 50 sweeps earlier. At around 



7 



70-120 sweeps there must exist a small instanton-like ob- 
ject that has been removed by the errors in the Symanzik 
action, but preserved by the tuned over-improved action. 

The second graph is a more typical example of what 
one sees when using the three different actions. The Wil- 
son action is still clearly the worst of the three, fluc- 
tuating the most. Meanwhile, the Symanzik and over- 
improved actions are fairly similar in their behaviour. 
Both stabilise at the same integer charge, but the over- 
improved action stabilises earlier. 

In the third graph we provide an example of how care 
must still be taken when using over-improved stout-link 
smearing. Contrary to the first two graphs, the curves in 
this graph seem to indicate that the over-improved action 
is the worst of the three. In this gauge field there must 
exist numerous instanton-like objects with size slightly 
greater than 1. Objects with this size will still be removed 
by over-improved smearing, but will survive for longer. 
Hence the topological charge takes longer to stabilise. 
The effect was great in this case because of the large 
lattice size, which meant it was possible for a few of these 
objects to exist on the lattice. The probability of finding 
such small objects on smaller lattices is significantly less. 

The final graph is a sample of a 16 3 x 32 lattice. It 
is shown here to demonstrate how it is much easier to 
smooth a smaller gauge field. Note that not all small 
lattices are this simple to smooth and we occasionally 
see behaviour similar to that in the top 3 graphs. In 
the larger lattices, the larger size means that there is 
a greater probability of finding an unstable topological 
object and it becomes more difficult to achieve integer 
charges. 



B. Topological Charge Density 

For the next part of the analysis we will directly ob- 
serve the topological charge density of the gauge fields. 
Our aim is to directly observe the differences in the gauge 
fields revealed by using the Wilson and over-improved ac- 
tions. 

To achieve this we will require a gauge field where the 
final topological charges from the two smearing proce- 
dures differ. We also consider a smaller 16 3 x 32 lattice 
because smaller lattices often provide clearer visualisa- 
tions. 

The topological charge, as a function of the number of 
smearing sweeps, is shown in Fig. [5] It appears as though 
an anti-instanton is being destroyed by the Wilson action 
from about 20 sweeps onwards. It will be interesting to 
visualise q{x) in this region to see if we can observe this 
behaviour. Indeed, by considering the differences in the 
charge density, we were able to locate the anti-instanton 
that is removed by the Wilson action. 

In Fig. [7] we show how the anti-instanton is affected 
by the Wilson action, and in Fig. [5] we have the corre- 
sponding charge density from the over-improved action. 
The pictures represent a single slice of the charge den- 



3 



-1 





1 1 1 1 
— e 


= 0.0 






— —0.25 




— € 


= +1 — 




i,i, 





10 20 30 40 50 

Number of Sweeps 



FIG. 6: The topological charge evolution under smearing for 
a 16 x 32 lattice. We see that when the Wilson action is used 
an anti-instanton is destroyed at around 20-40 sweeps. Some 
visualisations of the topological charge density in this region 
are provided later in the text. 



sity of the 4-D lattices as they evolve under the stout- 
link smearing. Objects with negative topological charge 
are coloured blue to green, and the positive objects are 
coloured red to yellow. 

After 30 sweeps we see that both smearing methods 
have revealed a similar vacuum structure. The effects 
of the errors in the Wilson action are first seen after 33 
sweeps, when the anti-instanton like object on the right 
begins to unwind in the upper-right corner. Here the 
charge density is approaching zero and therefore is not 
rendered. In a few sweeps the action density in this re- 
gion will manifest itself in the opposite winding, largely 
eliminating the total topological charge. The net effect 
is to suggest that the instanton-like object on the right 
invades the neighbouring negative object. However, the 
change in Q indicates that this is not an instanton - anti- 
instanton annihilation. At this point the majority of the 
negative topological charge density is lost and the to- 
tal Q for the configuration approaches 1. This kind of 
phenomenon should not be seen as filtering is applied to 
a lattice, and indeed it does not occur when using the 
over-improved smearing. 

After 36 sweeps the opposite winding has grown in 
size and it continues to grow in size as more smearing 
is applied to the lattice. After 39 sweeps the negatively 
charged object has all but disappeared. Although not 
shown, eventually the neighbouring positive object com- 
pletely engulfs the region originally occupied by the neg- 
atively charged excitation. 

This is a direct demonstration of how the discretisation 
errors in the Wilson action have resulted in an erroneous 
picture of the vacuum, and how by modifying these errors 
in the over-improved algorithm we are able to present a 
more accurate representation of the vacuum. 




FIG. 7: The evolution of the topological charge density for 
various sweeps of standard stout-link smearing. The sweeps 
shown are; 30 (top), 33, 36, 39 (bottom). In grey-scale (print) 
the darker areas represent regions of negative charge, and 
the lighter areas represent positive charge. In colour (online) 
blue to green represents negative topological charge and red 
to yellow represents positive. We see that a rather large anti- 
instanton is unstable under this smearing and is removed from 
the lattice, presenting an erroneous view of the vacuum. 
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FIG. 8: A visualisation of the topological charge density of 
the same gauge field shown in Fig. [7] this time with over- 
improved stout-link smearing. We see that in this case the 
anti-instanton in the lower right corner of the lattice is stable 
under smoothing, and remains stable for at least 200 sweeps. 
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VI. CONCLUSION 

We have demonstrated how to define an over-improved 
stout-link smearing algorithm, with the aim of preserv- 
ing instanton-like objects on the lattice. Using the new 
definition we showed how to select a suitable value of the 
parameter e, and suggest a value of —0.25. With the 
procedure defined, we demonstrated the success of the 
stout-link algorithm in preserving topological structures 
which were destroyed when using the standard Wilson 
action. The over-improved stout-link smearing can be 
used in future studies of vacuum structure or other simi- 
lar applications, where preserving topology on the lattice 



is important. 
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